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We consider collective laser cooling of atomic gas in the 
Festma lente regime, when the heating effects associated with 
photon reabsorptions are suppressed. We demonstrate that 
an appropriate sequence of laser pulses allows to condense a 
gas of trapped bosonic atoms into the ground level of the trap 
in the presence of collisions. Such condensation is robust and 
can be achieved in experimentally feasible traps. We extend 
significantly the validity of our previous numerical studies, 
and present new analytic results concerning condensation in 
the limit of rapid thermalization. We discuss in detail neces- 
sary conditions to realize all optical condensate in weak con- 
densation regime and beyond. 
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I. INTRODUCTION 

Laser cooling has led to spectacular results in recent 
years 0. So far, however, it has not allowed to reach 
temperatures for which quantum statistics become im- 
portant. In particular, evaporative cooling is used to 
obtain Bose-Einstein condensation of trapped gases . 
Nevertheless, several groups are pursuing the challenging 
goal of condensation via all-optical means . The in- 
terest in laser-induced condensation stems from the fact 
that, (a) a laser-cooled system is open and is not neces- 
sarilly in thermal equilibrium, i.e. its physics is in princi- 
ple much richer than that of evaporatively-cooled gases; 
(b) laser cooling is perhaps the only method which may 
condense all atoms into the ground state, i.e. reach an 
effective zero temperature. 

In traps of size larger than the inverse wavevector, 
fc^ 1 , the temperatures required for condensation are typ- 
ically below the photon recoil energy, Er = Hujji = 
Ti 2 k 2 L /2m, where m is the atomic mass. There exist sev- 
eral laser cooling schemes to reach such temperatures 
citeVSCPT, Raman. They exploit single atom "dark 
states" , i.e. states which cannot be excited by the laser, 
but can be populated via spontaneous emission. The 
main difficulty in applying dark state cooling for dense 
gases is caused by light reabsorption. Indeed, these states 
are not dark with respect to the photons spontaneously 
emitted by other atoms. Thus, at sufficiently high densi- 
ties, dark state cooling may cease to work, since multiple 
reabsorptions can increase the system energy by several 
recoil energies per atom [|-12|. In particular, laser in- 
duced condensation is feasible only if the reabsorption 
probability is smaller than the inverse of the number 



of energy levels accessible via spontaneous emission pro- 
cesses |10|| . 

Several remedies to the reabsorption problem have 
been proposed. First, the role of reabsorptions increases 
with the dimensionality. If the reabsorption cross sec- 
tion for trapped atoms is the same as in free space, i.e. 
~ 1/fcf,, the reabsorptions should not cause any prob- 
lem in ID, have to be carefully considered in 2D, and 
forbid condensation in 3D. Working with quasi-lD or - 
2D systems is thus a possible way to reduce the role of 
reabsorptions []l3|| . The most promising remedy against 
this problem employs the dependence of the reabsorption 
probability for trapped atoms on the fluorescence rate 7, 
which can be easily adjusted in dark state cooling |l4| . 
In particular, in the so called Festina lente limit, when 7 
is much smaller than the trap frequency u> |i"5| ], the re- 
absorption processes in which the atoms change energy 
and undergo heating are suppressed. This regime can be 
achieved either in Raman cooling Jl6| by changing the in- 
tensity and/or the detuning of the repumping laser, or in 
systems which simply posses a narrow natural linewidth. 

In a series of papers J!?] we have investigated col- 
lective cooling schemes in traps of realistic size in the 
Festina lente limit, and we have shown that laser in- 
duced condensation is possible. In this paper we present 
also such investigation taking into account the effects of 
atom-atom collisions. This paper is, in a sense, a final 
paper in a series of papers about the Festina lente regime. 
In the first two of those papers we have generalized dark 
state (Raman) cooling method beyond the Lamb-Dicke 
regime, i.e. beyond the regime when the size of the trap is 
smaller than the laser wavelength |L7]jlJ]. In particular, 
in the second of these references we have demonstrated 
the possibility of cooling to an arbitrary state of the trap. 
We have then applied those methods to samples of atoms 
and demonstrated possibility of all optical condensation 
first for the case of an ideal gas jl9) , and then taking into 
account elastic interatomic collisions [Q. In all of the 
above paper we have limited, however, our investigations 
to relatively small traps with Lamb-Dicke parameter r\ 
not greater than 5 in ID, and not greater than 2 in 3D 
because of numerical complexity. Also, we have consid- 
ered so far only the regime of weak condensation, when 
the energy levels of the trapped gas are not yet modi- 
fied by the mean field interactions. That implied that we 
could have considered only moderate numbers of atoms 
in the trap of the order of few hundred up to one thou- 
sand. 

In this paper we overcome the above mentioned short- 
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comings. In particular: 



II. QUANTUM MASTER EQUATION 



• Wc extend our numerical studies to traps with rj = 
8 in ID and n = 4 in 3D. 

• We present analytic description of the dynamics in 
the weak condensation regime and in the limit of 
rapid thermalization due to elastic collisions. 

• We extend our theory beyond the weak conden- 
sation regime to the hydrodynamic regime, and 
present the analytic description of the dynamics in 
the limit of rapid thermalization using Bogoliubov- 
Hartree-Fock theory at finite temperature. 

• Finally, we present a detailed discussion of require- 
ments that have to be fulfilled in order to achieve 
experimentally all optical condensation in the Fes- 
Una lente regime. 

The paper is organized as follows. The sections II-IV 
deal with the weak condensation regime. In section II, 
we formulate the Master Equation (ME) that describes 
the combined effects of dark state (Raman) cooling in the 
Festina Lente regime, and elastic atom-atom collisions. 
Cooling is dynamical, and consists of sequences of pairs of 
laser pulses inducing stimulated and spontaneous Raman 
transitions between the two electronic levels of trapped 
atoms, \g) and |e). The stimulated Raman transition 
induces the energy selective transition \g) — ► |e) that 
depopulates all motional states except the ground state 
which is dark; the spontaneous Raman transition via a 
third auxiliary level, similarly as standard spontaneous 
emission, is non-selective, and repumps the atoms from 
|e) to \g) populating all accessible motional states. The 
collisions introduce a thermalization mechanism which 
redistribute the populations of the trap levels according 
to a Bose-Einstein distribution (BED) . In section III we 
present our numerical results. We simulate the dynamics 
generated by the ME in ID and 3D, and show that laser 
induced condensation into the ground state of the trap 
is not only possible with collisions, but it is even more 
robust when they are present. In Section IV we present 
analytic results concerning the limit of rapid thermaliza- 
tion, in which the state of the system at each instant can 
be described by the thermal density matrix of the canon- 
ical ensemble. We demonstrate here that laser cooling 
leads to final temperature of the order of Er or less. 

In section V we extend our theory beyond the weak 
condensation regime. This can be done in the rapid ther- 
malization limit, since in that case one can describe the 
state of the system at each instant to be thermal and well 
described by the Bogoliubov-Hartree-Fock theory at fi- 
nite temperature T |2^| . Laser cooling introduces then a 
slow process of decrease of the effective temperature. We 
show that cooling down to T ~ is also possible beyond 
the weak condensation regime. This allows us to present 
a detailed discussion of the prospects for all optical BEC 
in Section VI. 



We consider N bosonic atoms with two levels \g) 
and |e) in a non- isotropic dipole trap with the frequen- 



cies UJ 



different for the ground and the ex- 
cited states, and non-commensurable one with another. 
This assumption simplifies enormously the dynamics of 
the spontaneous emission processes in the Festina Lente 
limit. We assume weak absorption pulses, so that no 
significant excited-state population is present. This al- 
lows to adiabatically eliminate the excited-state contri- 
bution, and consequently to consider the density matrix 
p(t) describing all atoms in the ground state \g), and di- 
agonal in the Fock representation corresponding to the 
bare trap levels. The ground state \g) into which the 
condensation will take place, should be the ground elec- 
tronic state of the atom, in such a way that the inelastic 
collisions of two \g) atoms are not possible ]22] ]. The Ra- 
man lasers are red-detuned sufficiently far from molecu- 
lar resonances, such that photoassociation losses can be 
neglected for the re gim e of atomic densities considered, 
n < 10 15 atoms/cm 3 [23]. For these densities, three-body 
losses are also negligible. We shall show that in the 
weak condensation regime, when the mean-field energy 
Ej < hui, the laser-cooling and collisional effects can be 
considered separately. In particular, the latter ones can 
be described by a quantum Boltzmann ME |25|,^6|]. The 
extension of our results beyond the weak condensation 
regime {Ej > Tiuj) is discussed in section V. 

In the following we follow the same notation as in Refs. 
p9| , p0[ . Let g m , g^ (e/, e\) be the annihilation and cre- 
ation operators of atoms in the ground (excited) state 
and in the trap level m (I), which fulfill the bosonic com- 
mutation relations. Using the standard theory of quan- 
tum stochastic processes one can develop the quantum 
ME which describes the atom dynamics in the Festina 
Lente regime 19 
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Here, 27 is the single-atom spontaneous emission rate, 
51 is the Rabi frequency associated with the atom tran- 
sition and the laser field, r/i m (k L ) — (e,l\e' lkL ' r \g,m) are 
the Franck-Condon factors, W(9, <p) is the fluorescence 
dipole pattern, uj® n (uf ) are the energies of the ground 
(excited) harmonic trap level m (I), and 6 is the laser 
detuning from the atomic transition |^7| . In the regime 
we want to study, only s-wave scattering is important, 
and then: 
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where ip m . denotes the harmonic oscillator wavefunc- 
tions, and a is the s-wave scattering length. 

In the following we are going first to work in the weak- 
condensation regime, where no mean-field effects are con- 
sidered. In typical experiments this condition requires 
that the condensate cannot contain more than 1000 par- 
ticles 26 1 . We shall also consider that <C oj. Due to the 
Festina-Lente requirements, Q is in general smaller than 
the typical collisional energy. Therefore we can formally 
establish the hierarchy £ ^> £ 2 > £\- 

Let us project into the ground state configurations 
with Nj atoms in the j-th level, and no excited atoms, 
\n) = \No,Nx, . . . ;<7)®|0, 0, . . .; e), using the projector 
operator VX = J^n \n) (n\(n\X\n) . Using Born-Markov 
approximation as in Ref. fl26|] , and expanding up to order 
0(£f), the ME for the dynamics of the reduced density 
operator v = V p becomes: 



v(t) = C coU v(t) +C coo iv(t), 



(5) 



where £ co u — — V £ 2 £§£ 2 describes the collisional part, 
and has the form of a QBME as that of Refs. H& 
The laser-cooling dynamics is described by 



C-cooi — —V£i[£o] 1 £\ 



dr / e -£o(r-r') £ie -£or' £i! (g) 



which has the same form of the ME calculated for 
ideal-gas case |lS| ]. The first correction to such split- 
ting between both dynamics is of the order £ 2 £\[£q]~ 2 , 
and therefore the independence between the collisions 
and laser cooling dynamics is only valid in the weak- 
condensation regime. In the limit of fast thermalization 
beyond the weak condensation regime, we will employ 
the hierarchy £ + £ 2 ^> £\- 



III. NUMERICAL RESULTS 

Eq. (|5|) can be simulated using standard Monte Carlo 
procedures. We have first performed such simulations for 
a relatively small sample of N = 133 atoms, and therefore 
well in the weak-condensation regime. Franck-Condon 
factors and trap frequencies can be efficiently approxi- 
mated using the states of an isotropic trap of frequency 
uj |Q. We concentrate on 3D cooling, whose analysis 
is greatly simplified by considering ergodic approxima- 
tion [ffil , i.e. considering that the population of degen- 
erate energy levels equalize on a time scale much faster 
than the collisions between levels of different energies, 
and than the laser-cooling time. The harmonic trap is 
considered beyond the Lamb-Dickc limit, i.e. the Lamb- 
Dicke parameter 77 = y/ En / Uui is larger than one. Due 
to memory storage and calculation times, our numerical 
simulation has to be restricted to not too large values of 
j]. In that case the atom (having initially an energy of 
the order of few Er) may increase its trap energy level in 
the spontaneous emission process by energies ~ En, and 
non-standard cooling schemes have to be used to avoid 
such heating effects jl7| . Generalization of the approach 
of Ref. (17) allows to cool dynamically (i.e. by changing 
absorption laser pulses app ropriately) individual atoms 
to arbitrary trap levels [|18| . 

Let us first consider the ideal gas case, i.e. a = 0. 
This approximation applies in principle only for very di- 
lute atomic samples. However, it has been demonstrated 
recently that a can be externally modified by magnetic 
fields Theoretical proposals describing modifica- 

tions of a using red-detuned lasers |24], or dc-electric 
fields [^9| are also under experimental studies. These 
methods could lead to experimentally feasible situations 
in which a ~ 0. In fact, very recent results of the JILA 
group indicate the realization of an "ideal" gas of 85 Rb 
is indeed possible |30(| . In such a case, it has been shown 
p9| that an appropriate dynamical cooling scheme would 
allow to condense a collection of trapped bosons, not only 
into the ground state of the trap, but also into excited 
ones. 

The cooling process that we consider consists of pulses 
of different frequencies, which are specially designed to 
leave just one single state not emptied. Since the level 
can be filled via spontaneous emission, the level acts as a 
trapping state, where the condensation is induced. Par- 
ticularly important is that it is possible to design two 
different dark-state mechanisms based, respectively, on 
the properties of the Frank-Condon factors, and on the 
destructive interference between the absorption ampli- 
tudes of orthogonal lasers ("interference" -dark-states). 
For a collection of atoms, the bosonic enhancement al- 
lows for faster cooling, the dark-states become more ro- 
bust with respect to changes of the physical parameters, 
and the multistability phenomena may appear. 

Using Monte Carlo methods, we have first analyzed 
the case of rather small trap, and studied the ground- 
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state cooling of 133 atoms in a 3D isotropic harmonic 
trap with 77 = 2, using 20 3D-energy shells (i.e. 1771 trap 
levels). We consider the case of 7 = 0.04a; and — 0.03w 
to assure the Festina lente conditions. This choice of 
parameters is rather restrictive. In practice, however, 
much larger values of 7 and fi, of the order of, say, u>/2 
can be used; in fact the use of larger 7 and Q is highly 
recommended in order to "hurry up slowly" , i.e. avoid 
long cooling times and reabsorptions simultaneously. 
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FIG. 1. Average of the population of the different energy 
shells, after evolving the system under collisions, starting from 
a thermal distribution of (n) = 6. 



In order to compare our results with the case of fi- 
nite a > 0, we employ as an initial atomic distribution 
the Bose-Einstein distribution (BED) (see Fig. 0) ob- 
tained after evolving the system in presence of collisions 
from an initial thermal distribution of mean 6hio. As 
we see, it already contains quite substantial amount of 
atoms condensed in the ground state, but also a lot of 
uncondensed ones. Wc must note that we begin with an 
BED below the condensation temperature due to numer- 
ical limitations, but qualitatively the same results could 
be found if the initial BED would not have been con- 
densed. We shall demonstrate this point later. The BED 
does not coincide exactly with the thermodynamical one 
described by grand canononical ensemble, due to finitc- 
size effects. Occupation numbers N n correspond to the 
ones obtained from the microcanonical ensemble at finite 
N and fixed energy. We apply our laser cooling cycles, 
each composed by two laser pulses of detuning S = sw, 
with S1.2 = -4,0, and time duration T = (2j)/Q 2 . The 
laser pulses are emitted in three orthogonal directions x, 
y and z, and are characterized by their respective Rabi 
frequencies Q x = Q y = O, Cl z = A z Vt. For the first pulse 
we assume A z = 1, while for the second one A z = — 2 is 
considered. With this choice, the second pulse is an "in- 
terference" dark-state pulse for the ground-state of the 
trap. Fig. |2|(a) shows (dashed line) that these two pulses 
are able to condense the population into the ground state 
of the trap, in absence of collisions; in particular no 
confinement pulses ]l7| , |l9| , ^8t (of detuning 8 = — 3r/ 2 w) 



are needed. This is due to the bosonic enhancement, 
and the fact that initially the system is already partially 
condensed. The dark-state pulse repumps the popula- 
tion in those states which are dark with respect to the 
pulses with detuning S = —Auj. During the dynamics the 
ground-state is the only trap level which remains not 
emptied. As pointed out in Rcf. the many body 

effects introduce one very important element to the dy- 
namics: Bose enhancement factors, that speed up the 
dynamics enormously. 
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FIG. 2. Population dynamics of the ground state of the 
trap for the case of 77 = 2, 7 = 0.04w, Q = 0.03w, TV = 133, 
starting from the BED distribution of Fig. [l]. (a) Cooling 
cycles composed of two laser pulses of detunings 5 — sui, with 
si,2 = —4,0, and A z — 1,-2 respectively, are used. The 
non-ideal (solid line) and ideal (dashed line) gas cases are 
compared. For uj — (2n) x 1kHz the cycle duration is ~ 
0.028s. This time can be, however, significantly optimized 
as discussed in the text, (b) Same case, but only employing 
pulse 1. 



Let us consider now the case in which collisions are 
taken into account. We shall consider the same situation 
and laser-cooling scheme as previously. Fig. ||(a) shows 
(solid line) the dynamics of the population of the ground- 
state in presence of collisions. After 600 cycles, all the 
population is transferred to the ground state of the trap. 
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This means that applying the laser cooling scheme brings 
the system into an effective BED of T = 0. It is easy to 
understand why the effect is maintained in presence of 
collisions, even considering that the collisional dynam- 
ics is much faster than laser-cooling. The laser-cooling 
mechanism tends to decrease the energy per particle (i.e. 
the chemical potential of the system), in the same way 
as evaporative cooling does, but without the losses of 
particles in the trap during the process. Thermalization 
via collisions brings the system to a lower temperature. 
If one repeats the laser cooling sufficiently many times, 
the system ends with an effective zero temperature. Fi- 
nally, let us point out that some auxiliary pulses which 
are needed in the ideal gas case, are not necessary in pres- 
ence of collisions. In particular, for the previous exam- 
ple, the pulse of zero detuning (required for the ideal gas 
case, Fig. ||(b) dashed line) is no more needed, as shown 
in Fig. ||(b) (solid line). Thus, efficient laser-cooling is 
not only possible in presence of collisions, but it can even 
be significantly simplified. 

We have been limited to rather small rj due to com- 
putational complexity. In particular, in simulations dis- 
cussed above 77 = 2 in 3D is considered. The maximal 
value of rj we could reach in 3D recently was 4, and the 
corresponding results are presented in Fig. 3. 
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FIG. 3. Population dynamics of the ground state of the 
trap for the case of rj = 4, 7 = 0.04w, Q = 0.03w, N = 500, 
starting from the BED distribution with mean energy 12 
Tilu. Cooling cycles composed of two laser pulses of detunings 
8 — slu, with Si j2 = —16,0, and A z — 1,-2 respectively, are 
used. The results for noninteracting (dotted line) and inter- 
acting (solid line) bosons are presented. For u> = (2n) x 1kHz 
the cycle duration is ~ 0.028s. This time can be, however, 
significantly optimized as discussed in the text. 

The Monte-Carlo simulations for 77 = 4 have been per- 
formed for 500 atoms occupying 33 3D-energy shells (i. 
e. 6545 energy levels). The initial mean energy per atom 
is equal to 12 hu>, which is above the condensation point 
(mean energy corresponding to T c is about 2.7tku). The 
initial number of condensed atoms is, however, relatively 
large (about 10%), what makes the cooling dynamics suf- 



ficiently fast from the beginning. The sequence of cooling 
pulses is similar as in the previous 3D calculations. The 
first pulse with detuning S — —16uj and A z = 1 is re- 
sponsible for confining the atoms, the second one with 
S = and A z = —2 is the dark-state pulse. This combi- 
nation of pulses is probably not the most efficient in the 
case of ideal gas, what may be noticed from Fig. [|. Af- 
ter about 200 cycles, the populations of the first excited 
levels become significantly large. Those levels begin to 
compete with the ground state and the cooling efficiency 
is reduced. This effect does not occur in the presence of 
interactions. The collisions distribute the atoms in more 
uniform way among the excited states, and the conden- 
sation is approached quite fast. 

2D and ID simulations show that same results hold for 
larger 77. In particular, in ID we have been able to study 
the cases up to 77 = 8, which is illustrated in Fig. f|. For 
the presented case, the values of parameters are chosen 
close to the limits of the applicability of our approach. 
This choice of the parameters leads to shorter cooling 
times, keeping Festina Lente condition. The pulse length 
is the same as in the previous simulations, in order to 
simplify the comparison of the cooling times. The pres- 
ence of interactions accelerate the cooling process, sim- 
ilarly to the considered 3D examples. The studies for 
the different values of Lamb-Dicke parameter allowed us 
to observe that our results scale with r\. In 3D, for in- 
stance, keeping N fixed the density scales as i]~ 3 , the 
three-body loss rate as ?7~ 6 , w as rf 2 , and the cooling 
time as rj 2 times the number of necessary cycles. The 
estimation of the latter is discussed in Sec. VI. 




cycles 

FIG. 4. Population dynamics of the ground state of the ID 
trap for the case of 77 — 8, 7 = 0.5u, £1 — 0.4a;, N = 500, start- 
ing from the BED distribution with mean energy 50hui. Cool- 
ing cycles composed of four laser pulses of detunings 5 = su), 
with si,2,3,4 = —64,-65,37,39 respectively, are used. The 
interacting case (solid line) and noninteracting case (dotted 
line) are compared. For ui = (2tt) x 1kHz the cycle duration is 
~ 0.056s. This time can be, however, significantly optimized 
as discussed in the text 
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IV. RAPID THERMALIZATION IN THE WEAK 
CONDENSATION REGIME 

From the above discussion it is clear that physically it 
is to be expected that collisions will act on a faster time 
scale then the laser cooling. In that case, collisions will 
provide a fast thermalization mechanism. The dynamics 
of the system will reduce to a fast approach toward in- 
stantaneous thermal equilibrium, and a slow change of ef- 
fective temperature T. In order to describe these changes 
quantitatively, we have to write down explicitly kinetic 
equations describing populations of the trap levels. This 
is done, as mentioned above, by performing adiabatic 
elimination of the excited state, and is discussed in more 
details below. 

In the regime of parameters we consider, the most of 
the atoms are in the ground state during the dynam- 
ics. In particular the system is affected by two different 
time scales, a slow one given by the commutator with 
Hi as and other given by the rest of terms in the master 
equation. Due to this fact we can use adiabatic elimina- 
tion techniques to remove the excited-state populations. 
Let \n) = \N 0} Ni, . . . ; g)<8>|0, 0, . . . ; e) the ground state 
configurations with Nj atoms in the jth level, and no ex- 
cited atoms; let their corresponding energies are Eft, and 
their corresponding populations pnn = (n\p\n). By using 
standard Projection Operator techniques it is possible to 
show that the adiabatic elimination of the excited state 
levels (see Appendix B of the Ref. p9|) leads to a set of 
rate equations for the populations: 

Pm(t) = Pn^,nPrnm.{t) ~ ^ Pffi*-nPnn{t) + f>m(t), 

rh rh 

(7) 

where P^ w -j are also defined in the Appendix of the 
Ref. ]l9j , and can be understood as the probabili- 
ties to undergo a transition from a ground-state con- 
figuration: |n) = \Nq, Nx, ■ ■ •) to other configuration 

|m)H^,iV{,...). 

The last term in Eq. (0) denotes the collisional con- 
tribution of the changes of Pnn(t)- A detailed discussion 
of the dynamics introduced by this term can be found in 
p5|p6[ . In our case, we do not have to specify it here: we 
just have to remember that in the regime of rapid ther- 
malization, this terms leads to a very fast thermalization 
of the system. 

Eq. O) can be rewritten in terms of rate equations 
with respect to the populations Nj in each level of the 
ground-state trap: 

N n =J2 T n^ m N m - T m ^ n N n + N™ 11 , (8) 

m m 

where the last term on RHS describes the contribution 
of collisions, while the rates are of the form: 

27 Jo Jo 
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[S - uj{1 - to)] + i"/R m i 



(N n + l-6 n , m ). (9) 



We have approximated here the excited and ground 
potentials by an harmonic potential of frequency w, 
whereas: 



flit pit 

Rmi = I dcj> d9W{9, 1 











Yhn>{k)\ 2 {N n ,+l-5 n ,, m ). 



(10) 



To obtain the rates T n ^ m we have used the properties of 
the creation and annihilation operators when applied on 
Fock states. Note that Eq. (|^) is the same as that we have 
found in |Q for the single-atom case. However, the rates 
@ and the single-atom ones are quite different. For the 
single-atom case they coincide as they should, but for the 
many-atom case, transition rates ^) become nonlinear, 
due to their dependence on the number of atoms in each 
trap level. In particular, we clearly distinguish here the 
two different quantum-statistical contributions: 

• In the denominator of the rates, we observe that the 
spontaneous emission acquires a collective charac- 
ter (similar as that observed in superradiance). 

• In the numerator of the rates, the bosonic-enhan- 
cement factor (N n + 1 — <5 n ,m) appears. 

The bosonic-enhancement factor favors the atoms into 
just one level of the trap. The reason is that if we are 
able to pump a significant amount of atoms into one sin- 
gle level \n), subsequent transitions into \n) are more and 
more probable. Therefore the condensation is more ef- 
fective. 

The contribution of the collective spontaneous emis- 
sion is unfortunately not so advantageous. The cool- 
ing methods designed for the single-atom case, which we 
want to apply in the many-atom case, are based on res- 
onant processes, and are therefore quite dependent on a 
narrow resonance. Note that the width of the resonance 
centered at 5 = u)(l — m) is given by ^R m i- If 7-RmZ 
grows, then the resonance is broadened, and the dark- 
state effects could cease to exist. On the other hand, as 
R m i grows, the height of the resonance becomes lower, 
and consequently the cooling becomes slower. These neg- 
ative effects can be avoided by using a sufficiently small 
7, which is possible in Raman cooling scheme. In Ref. 
|19[ for instance, the value 7 = 0.005a-> was used in all 
the calculations, which guaranteed for one-dimensional 
calculations that only the resonant terms were relevant. 
It has also been demonstrated that the negative effects 
are less important for higher dimensions. Note that the 
decreasing of 7 makes the cooling slower. In examples of 
R m is discussed in it has led to large cooling times 
(larger than 20 seconds for sodium atoms). This techni- 
cal problem can be eventually solved by decreasing ex- 
ternally 7 during the cooling process, in such a way that 



G 



the maximum effective coilective spontaneous emission 
rate remains approximately constant and comparable to, 
say, lo/2. Such technical optimization would increase the 
cooling rate, allowing for realistic cooling times, as we 
will show in Section IV. 

In the limit of fast thermalization, we can assume that 
all A n adjust their values to an instantaneous thermal 
equilibrium characterized by the temperature T(t), i.e. 
are given by the Bose-Einstein distribution. The tem- 
perature varies slowly in time as the laser cooling pro- 
ceeds. Particularly interesting is the situation below crit- 
ical temperature, in which we have 

N n = exp(-/?(t)7iw n )/(l - exp(-P(t)hu n )), (11) 



for n j^z 0, and 



A = A(l - (T(t)/T c ) 3 ), 



where (3(t) — l/ksT(t). The critical temperature is given 
by 

k B T c = hu(N/ 93 (i)) 1/3 , 

where the Bose function 53(1) ~ 1.2. From the rate 
equations (||) we immediately get the equation charac- 
terizing the slow changes of the bare energy of the sys- 
tem E = ^2 n hu; n N n . Below the condensation point, 
E{T) = 3fc B T(fc B T/?iw) 3 5 4(l), where g 4 (l) ~ 1.08, so 
that the above Eqs. || gives us the desired equation de- 
scribing the changes of the temperature, 



E = } TlUJ n 



^ ^ — m A m ^ Tm< — n A n 



We can easily rewrite this equation as 
dT 



dt 



F(T), 



(12) 



(13) 



where 



F(T) = 



T 



AE(T) 
^ hu! n 



( f?n — mA m ^ — nN n 



(14) 



Assuming that the ground state is perfectly dark, i.e. 
r m ^n = for all m, it is easy to see that this equation 
has a (stable) stationary solution T = 0. If the ground 
state is not perfectly dark, one obtains a finite stationary 
value T st . The stationary temperature is plotted against 
N in Fig. ||. If A is sufficiently large T st is very small, as 
we expect. The cooling rate during approach to T ~ is 
of the order of typical T$^ m , and becomes significantly 
slower only at very low temperatures for very large iV 
(compare the discussion in fll9|| ). 




20000 40000 60000 80000 1e+05 
N 

FIG. 5. Stationary value of the temperature T s t/T c as a 
function of N for the case of 77 = 4, 7 = 0.04oj, SI = 0.03a;. 
The same results in logarithmic scale are presented in the 
inset. 



The value of T st can be estimated assuming that the 
final distribution of N n will more or less balance the laser 
cooling dynamics. Similarly as in Ref. we can con- 
sider the limit when Nq ~ A, N m <C No for m ^ in a 
self-consistent way. In that limit the expected stationary 
values of N m are 



N„ 



-oAVrv 



We have r„«_ ~ 0(1), r ^„ ~ O(N ), so that N n ~ 
0(7/0;) in the stationary state. In this limit, the dy- 
namics consists of independent jumps — > n — > p9[ , 
and only the band of n's corresponding to roughly two 
recoil energies is relevant (n < 2Er/Tiuj = 2rj 2 ). From 
the expression 

m<2ri 2 

N m = N-N , 

m^0 



we estimate then 
Tst 



'm<2r] 2 

N m/N 

rn^O 



1/3 



which implies that k B T st ~ k B T c r] 2 ('y/u;N) 1 / 3 ~ 
(7/cj) 1/3 £; fl . As we see as A — * 00, T st /T c — > as 
I/A 1 / 3 . This estimation provides in fact only an up- 
per bound for T st , due to the fact that it does not take 
into account Franck-Condon coefficients. Numerical cal- 
culations for r\ = 4, 7 = 0.04a; and SI = 0.03a; presented 
in Fig. 5 clearly indicates then T st — > even faster as 
1/y/N for A — * 00. We shall show in the next section 
that similar conclusions hold beyond the weak condensa- 
tion. 
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V. BEYOND THE WEAK CONDENSATION 
REGIME 



the unitary Bogoliubov transform to the atomic quantum 
field operators 33 1: 



Let us first point out that beyond the weak condensa- 
tion regime, the mean-field energy cannot be neglected. 
This has two-fold consequences: (i) The levels of the trap 
are non-harmonic, i.e. they are not equally separated, 
because their energies become dependent on the occupa- 
tion numbers; (ii) the wavefunctions are different, and in 
particular the condensate wavefunction becomes broader 
(we consider here only a > 0). The fact that the energy 
levels are not harmonic any more, complicates the laser 
cooling, but the use of pulses with a variable frequency 
and band-width should produce the same results as those 
presented here. The point (ii) implies that the central 
density of the interacting gas is much lower than the one 
predicted for noninteracting particles, and therefore the 
dangerous regime of n ~ 5 x 10 14 atoms/cm 3 is reached 
for larger number of atoms than in the ideal case. For ex- 
ample, for the case of Magnesium, wavelength A = 600nm 
and 77 = 8, the mentioned density regime is reached for 
just TV = 3.5 x 10 3 , whereas in Thomas-Fermi approxima- 
tion the same is true for N = 6.6 x 10 5 . Below this num- 
ber, the interaction between the particles is dominated by 
the elastic two-body collisions considered in this paper. 
We show below that if our laser cooling scheme could be 
extended beyond the weak-condensation regime, laser- 
induced condensation of more than 10 6 atoms would be 
feasible. 

In general, to extend the quantum master kinetic the- 
ory beyond the weak condensation theory is a formidable 
task (compare [p5 26|). In general we do not know how 



the many body states of the trapped gas change in the 
course of the dynamics. The situation simplifies, how- 
ever, in the limit of fast thermalization. The collisions 
assure then that the instantaneous thermal equilibrium 
is achieved within a fast time scale. Such equilibrium can 
be very well described for a weakly interacting Bose gas 
using the Bogoliubov-Hartree-Fock (BHF) theory [ pi) . 



9k = I d 3 x 
9t 



~ f = / d 3 x 



u* k (x)^(x) + Vk{x)^(x) 
u k (x)&(x) + v k (x)^(x) 



(15) 
(16) 



where u k and v k fulfill the Bogoliubov-de Gennes equa- 
tions with the eigenenergy fiuik- In particular, go? So 
are associated with annihilation or creation of condensed 
particles. Using the bare states of the trap, the above 
relation can be rewritten as 



Sfc = ^2 [ u ln9n + v kn gl] , 

n 

9k = K«S™ + V *kn9n] , 



(17) 
(18) 



where Ufe„ = J d 3 x u k (x)ip n {x) ■ The above equations can 
be inverted to obtain 



9; 



Ukn9k ~ V kn gl 



k 

» = Y [ U kn9k - V 



kn9k 



(19) 
(20) 



The instantaneous equilibrium corresponds to a ther- 
mal state of the system described by the approximated 
Hamiltonian ^2 k huj k g k g k . In order to describe the cool- 
ing dynamics we have to express the quantum master 
equation of section II in terms of the BHF operators. 
We need to do it only for the terms responsible for laser 
cooling, since the collisional and free parts are responsible 
only for reaching the instantaneous equilibrium. Follow- 
ing the same steps as in previous section, we perform 
the adiabatic elimination of the excited state, carefully 
replacing particles by quasiparticles at each step. 



A. Bogoliubov-Hartree-Fock theory 

The BHF theory has various versions; for our purposes 
perhaps the best approach is the one valid for a fixed 
number of atoms developed by Gardiner ]3l| ] , and Castin 
and Dum Q. In those approaches one docs not break 
the phase symmetry, and avoid thus the time dependent 
effects associated with the spreading of the phase distri- 
bution ]33] |. Unfortunately, the number conserving ap- 
proach is technically tedious. It is therefore easier to use 
the standard BHF theory, with a broken phase symmetry 
and in the Popov approximation. In this case we have to 
neglect the effects of slow spreading of the phase distribu- 
tion. This approach will be used below. In the standard 
BHF theory at finite T, we introduce the quasi-particle 
annihilation and creation operators, which are related via 



B. Rate equations for populations 

Let |n) = \N ,Ni,...;g)®\0,Q,...',e) be the ground 
state configurations with N n quasiparticles in the nth 
level, and no excited atoms; let their corresponding ener- 
gies are and their corresponding populations pnn = 
(n\p\n). By using the similar Projection Operator tech- 
niques as in the previous section it is possible to show 
that the adiabatic elimination of the excited state levels 
leads to a set of rate equations for the populations: 

Pnn(^) ^ Pn< — mPmm(^) ^ Prh< — nPnni^) ~t~ Pfin (0> 
m m 

(21) 

where Pm^n have similar form as those that appear in 
the weak condensation limit 
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Pn^ 1% = 2 7 / dcf>d9 sineW (0, </>) 



n\(fjf s (k)gl + Ci s (k)g s )ei 



Erf, — H, 



eff 



(22) 



and 



where fjuQe) = J d s xipf(x)e ik ' s u s (x) = rji s >u s 

Ch(k) = J d 3 xtpl(x)e ik - s v s (x) = J2 S > Vls' v ss' ■ The laser 
Hamiltonian is also reexpressed as 



Hias = ? yY^ e l(vim(k L )g m + (? m (k L )g f m ). 



As before, the probabilities (]22j) can be understood as the 
probabilities to undergo a transition from a ground-state 
quasi-particle configuration: |n) = \Nq,Ni, . . .) to other 
quasi-particle configuration |m) = |7Vg, N[, . . .). 

The last term in Eq. (21) denotes the collisional con- 



tribution to the changes of pnait). We do not have to 
specify it here: as in Section IV we just have to remem- 
ber that in the regime considered, this terms leads to a 
very fast thermalization of the system. 

Eq. ( pl| ) can be rewritten in terms of rate equations 
with respect to the quasi-particle populations N n in each 
level of the ground-state trap: 



rN„ 



,.iV» + N. 



coll 



(23) 



where the last term on RHS describes the contribution 
of collisions, while the rates are of the form: 



r 

i n< — r 



02 

27 



d(j> / dewie,^) 



ifii n ( k )mm(k L ) 



^ s 

i 



llm{k)Q m {k L ) 



(N n + 1- 6 n>m ) 



^8 

i 



(24) 



In the above expression we have used rotating wave ap- 
proximation with respect to quasiparticle frequency spac- 
ings. One should stress that, although in general split- 
tings between the quasiparticle energies can be small, in 
the most interesting limit of large N, i.e. in the hydro- 
dynamic regime the quasiparticle spectra are very regu- 
lar, and the splittings are of the order of the bare har- 
monic oscillator frequency u> p4| . That is why in the 
hydrodynamic regime the use of RWA with respect to 
the quasiparticle energy splittings is very well grounded. 
The widths in Eq. (p4) become now 



2k 



d(f) / d8W(0,(f>) 



(25) 



We have additionally included in the expression (]24|) the 
width 7l which mimics the effects of the laser bandwidth, 
and which is necessary to assure approximate fulfilling 
of the resonance conditions. The widths 7l should in 
practice be of the order of u>. 

The equations derived in this section have a very sim- 
ilar form to the one derived in the previous sections. In 
particular, we expect that in the Thomas-Fermi (hydro- 
dynamic) regime, they will lead to the self-consistent 
equation for the effective temperature. We can easily 
rewrite this equation as 



dT 
dt 



= nn 



(26) 



where 



F(T) 



E'(T) 



- n N n 



, (27) 



where N m are distributed according to the Bose-Einstein 
distribution corresponding to the quasiparticle Hamil- 
tonian, whereas E(T) = ^„ftwJVn) E'{T) = dE/dT. 
Note that temperature dependence enters the above 
equations explicitly, and through the dependence of the 
quasiparticles and their energies on the temperature in 
the BHF theory at finite T. In the limit of large N, 
the same argumentation as in previous section can be 
applied. We expect that the system will be cooled 
down to the temperature T st ~ 0, and that the sta- 
tionary temperature T st will be of the order of fesX^ ~ 
kBTc^^/tuN) 1 / 3 ~ (j/uj) 1/3 E r . As we see as N -> oo, 
T st /T c -> at least as fast as l/N 1 / 3 . 



VI. PROSPECTS FOR ALL OPTICAL BEC 

According to our results the prospects for all optical 
BEC are very good provided several precautions are re- 
alized in experiments. The prescription to achieve all 
optical BEC is: 

• Use a dipole trap of the moderate size rj ~ 2 — 8 
and oj ~ 2ir x 1 kHz. 

• Trap and cool atoms into the electronic hyperfine 
ground state. This allows to eliminate non-elastic 
two body collisions from considerations. 

• Work in the Festina-Lente limit to avoid reabsorp- 
tions. Use either natural narrow line, or quenched 
(Raman) transition of width J e ff < w. In fact 
to shorten the cooling dynamics try to work with 
Jeff ^ u/2. 

• work with red detuned laser tuned below or in be- 
tween the molecular resonances to avoid photoas- 
sociation losses. 







• Avoid high densities, i. e. limit the number of 
trapped atoms N in such a way that n < 5 • 
10 14 /cm 3 . This will allow you to avoid 3-body 
losses and remaining photoassociation losses. 

In the Tab. | we present our estimates for the max- 
imal number of condensed atoms and the cooling time. 
We have considered here Mg atoms, and assumed Al = 
600nm and scattering length a sc ~ 5 nm. For various 
values of r\ we have then calculated corresponding trap 
frequency u>. Maximal number of atoms N max has been 
estimated from requirement that the corresponding peak 
density should be 5 • 10 14 /cm 3 . The calculations has 
been done for the ideal (a sc — 0) and weakly interact- 
ing gas, where for simplicity the Thomas- Fermi density 
profile was assumed. 

For the estimation of the cooling time t coo i we use 
the following reasoning. The probability of atom jump 
(transfer) from level m to level n is proportional to the 
occupation numbers: p m —tn ~ N m (N n + 1). If we omit 
the dependence on Frank-Condon coefficients, the jump 
probability is given by p m _> n = A/'iV m (A r „ + 1), where Af 
should be determined from the normalization condition: 
E m E n fa-»n = Ptot, where p to t is total probability of 
jump. Below we shall assume that in each cycle roughly 
10% of atoms undergo a jump, i. e. p to t —0.1. We con- 
sider the final stage where only the dark-state pulse for 
the ground state is responsible for cooling. We assume 
that the ground state is perfectly dark, i. e. po— >m = 
for all m ^ 0. The spontaneous emission process scat- 
ters the atoms within the energy band of width about 
Arfhu) = AEr. Assuming that the dynamics takes place 
within the region of trap levels with energies smaller than 
4r] 2 huj, the normalization coefficient reads 



N = 



Ptot 



(N — N )(N + N st ) 



(28) 



where N st is the number of states with energies between 
and Ari 2 huj: N st « (4r; 2 ) 3 /6. The cooling efficiency grows 
if we assume that in each cycle 10% of the excited states 
population is transferred. This number is consistent with 
the conditions for the adiabatic elimination used in our 
derivations. The population of the ground state is in- 
creased in every step by (N — N Q ) J2 m ^o Pm-^0- Thus we 
are able to write the formula describing the ground state 
dynamics: 

N (k + 1) = No(k) + e(N - N (k))(N (k) + 1), (29) 

where e = (10(A + 7V st )) _1 and k labels the cycles. In 
rough approximation the equation (^9|) may be replaced 
by the differential equation 



Wo 
dk 



= e(N-N )(N + l). 



(30) 



Starting with 1 atom in the ground state, the number 
of cycles needed to obtain N/2 condensed atoms may 
be calculated to be ln(N)/(eN). Finally assuming that 



each cycle have the duration 2/7, where 7 ~ we 
estimate t coo i ~ 81n(iV)/ (uieN). The analytic estimates 
agree quite well with numerical simulations of Sec. III. 
Direct inspection of the Tab. |, clearly shows that all- 
optical BEC of quite large number of atoms should be 
feasible within a reasonable cooling time. 
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a sc = 






5nm 


V 


max 


tcooi 


N max 


tcooi 


2 


55 


119ms 


161 


59.3ms 


4 


439 


5.45s 


1.0 x 10 4 


441ms 


6 


1.4 x 10 3 


51.86s 


1.2 x 10 5 


1.21s 


8 


3.5 x 10 3 


232.4s 


6.6 x 10 5 


2.5s 



TABLE I. Maximal number of condensed atoms and cool- 
ing times calculated for Mg atoms, assuming Al = 600nm 
and scattering length a sc ~ 5 nm 
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